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^r '■ Abstract 

r^ . 

'^ ' Fifth order, quasi-linear, non-constant separant evolution equations 

cn . are of the form ut = A^-^ + B, where A and B are functions of x, t, 



u and of the derivatives of u with respect to x up to order 4. We use 
the existence of a "formal symmetry" , hence the existence of "canonical 
conservation laws" /?(j), i = — 1,...,5 as an integrability test. We de- 
fine an evolution equation to be of the KdV-Type, if all odd numbered 
canonical conserved densities are nontrivial. We prove that fifth order, 
C^ ' quasi-linear, non-constant separant evolution equations of KdV type are 

polynomial in the function a = A^'^; a = {au^ + /3'U3 +7)~^ > where a, 
f3 and 7 are functions of x, t, u and of the derivatives of u with respect 
to X up to order 2. We determine the M2 dependency of a in terms of 
P = 4a7 — /3^ > and we give an explicit solution, showing that there 
are integrable fifth order non-polynomial evolution equations. 

Keywords: Evolution equations, Integrability, Classification, Recur- 
sion operator. Formal symmetry 

1 Introduction 

In the literature, the term "integrable equations" refers to nonlinear equations 
that can either be transformed to a linear equation or solved by an inverse 
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spectral transformation. The methods to determine whether a given equation 
is hkely to be integrable or not are called "integrability tests". An extensive 
review of the integrability tests and a comprehensive list of integrable equation 
have been given in [1]. 

The Korteweg-deVries (KdV), Sawada-Kotera and Kaup hierarchies are well 
known hierarchies of integrable equations; the KdV hierarchy appears at all odd 
orders, while the Sawada-Kotera and Kaup hierarchies exist at odd orders that 
are not multiples of 3. Most of the constant separant third order integrable 
equations listed in [1] are transformable to an equation belonging to one of these 
hierarchies; it is in fact conjectured that the only exception is the Krichever- 
Novikov equation [2]. The search for new hierarchies of integrable equations 
starting at order seven didn't also give any positive result [3]; this situation 
was clarified for the polynomial and scale invariant evolution equations by 
the remarkable result of Wang and Sanders [4], stating that in the class of 
polynomial scale invariant equations, any integrable equation of order 7 or 
higher is a symmetry of a lower order equation. Thus, the open problems 
were reduced to the classifications at the third and fifth orders and to the 
classification of higher order non-polynomial equations. 

In a series of papers we have applied the "formal symmetry" method of 
[1] to generic, non-polynomial equations; we have first proved that integrable 
evolution equations of order 7 and greater are quasi-linear [5]; then we have 
shown that they are polynomial in top three derivatives [6]. Motivated by these 
polynomiality results, we decided to investigate the structure of fifth order 
quasi-linear integrable equations with non-constant separant to see whether 
they would be polynomial and transformable to the constant separant case 
classified in [1]. 

In the present work we consider quasi-linear fifth order evolution equations 
Ut = A^ + B and we focus on the case where p(3) is nontrivial, a case naturally 
distinguished by the classification algorithm using canonical conserved densi- 
ties. Evolution equations with an unbroken sequence of canonical densities will 
be called of "Korteweg-deVries (KdV)-Type" as opposed to the Sawada-Kotera 
and Kaup equations [7] for which p(3) is trivial. 

We introduce our notation and give basic definitions in Section 2. The 
dependency on ^4 is easily determined as presented in Section 3.1 and it is 
shown that Ut is quadratic in U4. The form of the dependency on u^ depends 
on whether the canonical conserved density p(3) is trivial or not. In Section 3.2.1 
we show that for the KdV- Type equations, the non-triviality of p(3) implies that 
A-*^/^ = a = (aul + I3u^ + 7)"^/^, where a, (3 and 7 are certain functions of x, t, 
M, Ml and U2- We use the conserved density conditions to determine the form 
of the evolution equation; in particular, in Section 3.2.3, we show that they are 



polynomial in a and in M3, by proving that the coefficients of logarithmic and 
other transcendental functions of a vanish by virtue of the partial differential 
equations satisfied by a, /3 and 7. The discriminant of the quadratic expression 
au\ + /3-U3 + 7 is denoted by —P and assumed to be negative, to avoid real 
roots. The U2 dependency of the functions a, /3 and 7 are determined in terms 
of P but the Ml, -u, X and t dependencies could'nt be obtained. In Section 3.2.4, 
we present a fifth order integrable equation, assuming that ut has no x, t, m. 
Ml dependency, hence we present a new integrable equation in this class. 

As discussed in Section 4, these results show that at the fifth order there 
are non-constant separant integrable equations which are possibly higher order 
analogues of the essentially nonlinear class of equations at the third order [8]. 
Furthermore, as part of ongoing work, we have shown that there are candidates 
of integrable equations of the form Ut = a^'^^^U2k+i + • • • at orders /c = 3, 4, 5, 6, 
which supports the existence of a new hierarchy, which nevertheless is expected 
to be transformable to the KdV equation by a Miura type transformation 
involving third order derivatives. 

2 Preliminaries 

2.1 Notation 

We shall work with evolution equations in one space dimension x and use the 
following notation for the partial derivatives of the dependent variable u. 



du du d'^ 



u 



ot ox ^-w^ ox" 



k~~times 



In the following we shall assume that all functions are sufficiently differentiable. 
Infinitely differentiable functions of x, t, u and of the derivatives of u up to an 
arbitrary but finite order are called "differential functions" [9]. We shall use 
subscripts for partial derivatives with respect to the derivatives of m, as below 



^k- TT-^ Vkj 



duk dukduj 

The labels of functions will be denoted by indices in parenthesis, such as (T(j) 
or p(i). 

A scalar evolution equation of order m is of the form 

Ut = F{x,t,U,Ui,...,Um)- (1) 

The total derivative operators with respect to x and t are respectively denoted 
by D and Df. That is, if (y9 is a differentiable function of x, t, u and the partial 



derivatives of u with respect to x up to order k, then 

dip dip dip dip 
-UV' = -^ + ^ Ml + -^ U2 + ... + ^ — Mfc+i = 'Px + Vo'^i H \- VkUk+i- (2) 

If u{x,t) satisfies (1) and ip is as above, then 

Dt^=^ + ^F+^DF + ... + ^D'F = ^t + ^oF + --- + VkD'F. (3) 
dt du dui duk 

2.2 The "Level Grading" 

The notions of "level grading" and level homogeneity introduced in [6] are 
analogues of scaling and scale homogeneity for polynomials; in the level grading 
context, the "level above a base level A;" reflects the number of differentiations 
applied to a function that depends on the derivatives of order at most k. The 
crucial property of the level grading is its invariance under integrations by 
parts. This property implies that "top level" terms of a conserved density p 
will give top level terms in its time derivative Dtp., hence we may omit lower 
level terms and still track correctly higher level integrability conditions. The 
"level of Un above the base level /c" is defined as n — k and clearly it is equal 
to the number of differentiations above k. The level of a monomial u'!} . . . u^T 
is then defined as hi{ni — /c) + • • • + bm{nm — k). For example, if the base level 
is /c = 4, -Us has level 1, while if the base level is A; = 3, -Us and u\ are both 
monomials of level 2, hence, for example, the right hand side of (14) is a sum of 
level homogeneous expressions of levels 2, 1 and above the base level /c = 3. 
In the present work, our results are based on explicit computations including 
derivatives of all orders, however we still indicate level related results whenever 
necessary. 

2.3 Formal Symmetries 

The Frechet derivative of a differential function F is defined by 

m rs-p m 

j=0 j=0 

Let a = a{x, t, u, ..., «„) and p = p{x, t, u, ..., «„) be differential functions. If 

(Tt = F^a, (5) 

then a is called a "symmetry" . If one can find a differential function rj such 
that 

pt = Dt] (6) 



then p is called a "conserved density" . Let i? be a linear operator. If for any 
symmetry a, Ra is also a symmetry then R is called a "recursion operator" . If 
Ra is a symmetry, then the equation 

{Ra\ = F,{Ra) (7) 

gives {R(j)t = RtCf + RF^,a = F^Ra, hence 

{Rt + [R,F,])a = Q. (8) 

Therefore, if R satisfies the operator equation (8), it certainly satisfies (7). 
Note however that in this case there is no guarantee that the quantity Ra 
which satisfies the symmetry equation (7) is a local function. This suggests 
that if we would solve the operator equation (8), we could as well work with 
a pseudo-differential operator R. In addition, we may only be interested in 
solving the operator equation (8) up to a finite order. A pseudo-differential 
operator that solves (8) up to a finite order is called a "formal symmetry". 
The existence of a formal symmetry has been proposed as an integrability test 
in [1]. This approach has the advantage that the existence of the recursion 
operator, hence the classification of integrable equations is insensitive to the 
order and the functional form of the operator R. We have used the symbolic 
programming language REDUCE in order to compute the partial differential 
equations implied by these conserved density conditions. 

3 The form of integrable equations 

The criterion of integrability by the existence of a formal symmetry [1] implies 
the existence of certain conserved densities p(j), i = —1,0,1,.... It is well 
knowu that for auy order m. if F„. = £. F„._, = ^ the,r p,_.. = 

Fm and p(o) = Fm-i/Fm are always conserved densities. The expression 

of the conserved densities p(^i), i = 1, . . . , 5 are given in Appendix A[10]. 

3.1 Dependency on U4 

We assume that the evolution equation is of the form 

Ut = a^U5 + B, (9) 

where a and B are functions of x, t, u, Ui, i = 1,...,4. Substituting (9) 
in the conserved densities given in Appendix A, we can obtain explicitly the 
form of the canonical conserved densities p(i) and p(2) . However the conserved 
densities p(3), p(4) and p(5) involve the expressions o"(_i), (T(o) and (T(i) defined 



by DtP(i) = Da(i), hence they can be exphcitly determined only after the lower 
order conserved density conditions are solved. Note that, since cr(3), (7(4) and 
(T(5) are not used anywhere else, we may omit the total derivatives in p(j), for 
i = 3,4,5. In practice we solve the lower order conserved density conditions 
up to a certain order, replace the results in the other conserved densities and 
proceed iteratively. After substitutions and integrations by parts, the conserved 
densities are of the form given below. 

P(-i) = a^\ 
P(o) = ba'^a^u^ + B^, 

P(l) = P{l)ul + Q{l)U5 + R{1), 

P(2) = M^2)UI + P(2)UI + 0(2)^5 + ^(2) , 

P(3) = ^(3)^6 + Q{3)ut + R(3)ul + ^(3)^5 + T(3)M5 + f/(3) , 

P(4) = M(^i)ulu5 + P(i)ul + Q{i)ul + R(i)ul + S'(4)M5 + T(4)M5 + U(^i)U^ 

+ V(4), 
P(5) = P{b)u] + Q{^)ul + R(5)ulu5 + S(5)ul + T(5)M^ + U(5)ut + V(5)ul 

+ iy(5)M2+X(5)M5+y'(5). (10) 

In the equations above B^ = ^ and the coefficients P(i), Q(i), R{i) {i = 
1, ...,5), %, f(i), f7(i) {i = 3,4,5), M(2), M(4), %), ^"(5), W^(5), -^^(5) and 
Yj^s) are functions of x, t, u, ui, U2, u^, u^ that can be computed from the 
expressions of the canonical densities given in the Appendix A. We omit their 
explicit expressions except for 

P(i) = -, P(3) = - a[a44a - 404]. (11) 

In [5], Corollary 4.3, it has been shown that conserved densities of order n 
larger than the order of the evolution equation should be at most quadratic in 
their highest derivatives; thus in particular, we should have M(2) = 0. 
The coefficient of u^u^ in Dtpi gives 

10040^(0440 - 404) = 0. (12) 

Thus we should either have 04 = 0, or if 04 ^ 0, then 0440 — 404 = 0. If 04 = 0, 
then the coefficient of u\ in DtP(^-i) is an expression that involves o and -B4. 
Differentiating this expression with respect to 1x4 twice we obtain P444 = 0. 

On the other hand if 0440 — 4o| = 0, then, P(5) = — 20o|o^ and (5(5) = 
^a^a^ and the coefficient of ufuQU^ gives 04 = 0. When 04 = 0, in p(4), 
M(4) = P(4) = Q(4) = -R(4) = 0. By the Corollary 4.3 of [5], p(4) should be 



quadratic in u^, hence 5(4) should be zero. This condition gives again ,8444 = 0. 
Thus in either case, the M4 dependency is determined as below. 

04 = 0, ^444 = 0. (13) 

3.2 Dependency on M3 

As a result of (13) the evolution equation reduces to the form 

ut = a^U5 + Buj + Cu4 + G (14) 

where a, B, C and G are now functions of x, t, m, Mi, U2 and ^3. 

With this form of F the canonical conserved densities reduce to the form 
below. Although the explicit expressions of the conserved densities have been 
obtained, we use the generic expressions given below for computational effi- 
ciency. 

P(-i) = a"\ 
P(i) = P{i)u\ + Q(i)U4^ + Rii), 

P(3) = {P{2)ul + Qi2)ul)+ R{2)UI + S(2)u\ + T(2)Ua + U2, 

P(5) = (^(3)^6 + Q(3)W5 + %)M5M4 + S{^)ul) + (T(3)M^M4 + U('i)u\) 

+ (^(3)^4 + W^(3)M^) + X(3)m3 + Y(^)UI + Z(3)M4 + M(3) (15) 

The coefficients P(j), Q(i), R{i) {i = 1,2,3), %), T(i), f/(i) (i = 2,3), V(3), W^(3), 
X(3), Y(3), Z(3) and M(3) are functions of x, t, u, ui, U2, M3. 
Definition 3.1. An evolution equation ut = F{x,t,u, . . . ,Um) of odd order 
m is called of "KdV-Type", if it admits an infinite sequence of non-trivial 
conserved densities at all orders. 

In terms of canonical densities the condition above means that all the odd 
numbered canonical densities P(2k+i), in particular, the canonical conserved 
density p(3) is non-trivial. Note that these conserved densities are respectively 
of levels 0, 2, 4 and 6 above the base level 3, hence their total derivatives with 
respect to time will be respectively of levels 5, 7, 9 and 11. 

We shall use the conditions obtained by equating to zero the coefficients of 
top 3 levels terms in the conserved density conditions as listed in Table 1. The 
coefficients of these monomials denoted as Kijk, where the first index i denotes 
the label of the conserved density, the second index j denotes the level of that 
term (as a hexadecimal number) and k counts the number of such conditions. 
For example, K^u is the 4th condition obtained as a coefficient of a term of 
level 11 (6 as hexadecimal) in DtP{z). 



Table 1: The conserved density conditions for p(j) getting from the coefficients 
of terms of top 3 levels. 



P{i) Top level Top-1 level Top-2 level 

031 
052 "4 " -'^042 



P(_i) ulu^ - Ko5i, ul -Ko4i, ul - K( 

■u5 - Kn'^0 ui - Kn 



P(i) ului - Km, ul - Kifsi, 



UlU4^ - Kn2, Ms - KiQ2, 

ulul - Kn^, ulul - 7^163, M5M4 - ^151, 

ul - Kiu ul - KiQ4, ul - Ki52 



P(2) UjUi - K291, 


M? - K28I, 


- 


ul - K292, 


- 


- 


ulu^Ui - K293, 


Ulu5 - K282, 


ului - K271, 


ulul - K294, 


ulul - K283, 


- 


ui,U4 - K295, 


4 jT^ 

M5 " -"-284, 


- 


ulul - K29&, 


■U5-U4 - K2S5, 


M5M4 - -ft'272, 


ulul - K297, 


2 4 Ti^ 
U^U^ - A 286, 


ulul - K273, 


ul - K29S 


ul - K2S7 


ul - K274 


P(3) ului - K^bl, 


Ug - K^al, 


- 


UjUe - K3b2, 


UjU^ - Ksa2, 


- 


u'^U5U4 - K^bS, 


- 


- 


UjUl - KsbA, 


UjU^ - K^aS, 


^7^4 - -f^39l! 


Ulu5 - K^ibb, 


- 


- 


ulul - Ksb6, 


UIU4 - KsaA, 


ul - K392 


ulului - K3b7, 


ulul - Ksa5, 


- 


ulur^ul - K^bs, 


Ulu^ul - K^aG, 


UIU5U4 - i^393, 


ulul - K3b9, 


U^U^ - Ksa7, 


ulul - Ks94, 


U^Ua - Ksba, 


Ur, - K^aS, 


- 


ujul - Ksbb, 


U^U^ - Ksa9, 


U^U4 - K395, 


ulul - Ksbc, 


U^U^ - K^aa, 


ulul - K39Q, 


ulul - Ksbd, 


U^U^ - K^ab, 


ulul - Ks97, 


Ul^ - Ksbe, 


M4° - Ksac, 


ul - Ks98 



3.2.1 Top level computations 

We start by collecting those top level equations that will determine B and give 

a third order ordinary differential equation for the M3 derivative of a. Although 

the explicit forms of the P(j) 's are known, it will be more convenient to consider 

them as new unknown functions. 

From the condition i^osi we get 

4 8 

^051 : 03330^ - Taggosa^ - - as^aB + Sa^a"^ - - alB = 0. (16) 

Recall that as we are studying KdV-Type equations we are assuming that p(3) 
is nontrivial, hence P(2) 7^ 0. We then use the equations 

K.M : 35a3a%)-5(P(3))3a5 + 45P(3)=0, 



^3bl 

Km 
K291 
K292 



15a3a'P(i) - 5(P(i))3 a' + ABP^,) = 0, 
25a3a^P(2) - 5(P(2))3 a' + 4PP(2) = 0, 
20a3a^P(2) - 5(P(2))3 a' + 6PP(2) = 0. (17) 

Taking the difference of the last two equations, we solve B as 

B = ^ a^a\ (18)r 

From the remaining equations in (17) we find that the derivatives of P(j) and 
a with respect to u^ are proportional, and we solve them as below 

P{i) = P(io)a\ P(2) = P(2o)a', P(3) = P{^o)a'. (19) 

where P{^io) are functions of x, t, -u, Ui and U2- Finally, i^osi gives 

03330^ — 8033030 + I203 = 0. (20) 

Assuming a = Z^ and substituting in (20), we see that for k = 1/2, Z333 = 0, 
hence o is of the form 

a= (aM32 + /3M3 + 7)-| (21) 

where a, (3 and 7 are functions of x, t, u, Ui and U2- From (21) we can see that 

1 1 

03 = --o^(2aM3 + /3), Oj = --o^(aiM3 + /3iM3 + 7j), i = a;,0,l,2. (22) 

If we require o to be real for every ^3, we see that 

/3^ - 4q;7 < 0, and a > 

and we define 

P = 4a7 - /32. (23) 

From the expressions, -ft'293, K^h2^ K^hi, and -^'368, we solve respectively (5(2), 
(5(3), P(3) and S'(3) but we omit the explicit expressions here. The rest of the 
top level conditions are automatically satisfied. 
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3.2.2 Top-1 level computations 

By computing the coefficients oi u1, uf and -Ug respectively in p(i), p(3) and p(5), 
we find that 

P(lo) = ^P, P(2o) = -gP, P(3o) = ^P. (24) 

From the (Top-1) level condition K2S1, C is obtained as, 

5 DP 
C = --—a^ + 5{Da-a3Ui)a*. (25) 

We can again solve R{2), T(^3) and V(3) respectively from K282, -^3a2 and K^ae 
but we omit these expressions. 

Then from K1Q2, the coefficients of m| we obtain 

a = a(o)P', (26) 

where q;(o) = a(o){u,ui), and the coefficient of ^3 gives 

^ _ i B—- - 3 p (— —u —^ r27) 

(9^2 2 du2 P \dui du dx J 

The remain conditions of the (Top-1) level conditions are satisfied automati- 
cally. 

3.2.3 Top-2 level computations 

All of the (Top-2) level conditions except for -ft'292 are identically satisfied. From 
the condition ^^292 we get a second order ordinary differential equation for the 
M3 derivative of G. The form of this equation suggests that we have to multiply 
it by a~^ in order to integrate it. After substituting 

{2au3 + P) = -2a3a-^, (28) 

the resulting expression is 



a-3 



6-33 - 80-^03^3 + Y. (^(0^3 + m) «"' = 0' (29) 



j=i 



where the A(j)'s and p(i)'s depend on x, t, u, Ui and U2. In order to integrate 
this equation, we need to know the integrals of various powers of a. Note that 
a~^ is an irreducible quadratic polynomial in 1x3, hence the integrals of a^u^ and 
a? will involve logarithms and arctangent functions. However, the coefficients 
of these terms will be zero, hence, instead of writing these integrals explicitly. 
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we just indicate these as new functions. These integration formulas are given 
below. 

aus dus = / a du^, 

aa 2a; J 

f a^ du3 = i) = 2P-1/2 ^g^^-i (p-i/2(2aM3 + P)) , 

1 P [ 

a?u^ du3 = if = ln(a) a^ dus, 

a 2a J 

2 

a^ du3 = — {2au3 + /3) a. (30) 

The integrals of a'' and a'^u^ for A; > 3 are evaluated using the iterative 
formulas given below. 

«' du3 = T^]s (2«% + /3) a'-' + ^^^^^ a I a'-' du^, (31) 
K — 2 F k — 2 F J 



a^us du3 = — ( hi) a""^ - -— j a" du^. (32) 



The resulting integral is of the form 



3 



G^a ^ + /t(i)^ + /t(2)V9 + J^ (z^(i)M3 + V(i)) a^' = X, (33) 

where the fi;(i)'s, z/(j)'s, 77(j)'s and x depend on x, t, -u, ui and ^2- The coefficients 
of the functions ip and ip defined in Eqn.(33) are zero, after we substitute the 
equations for a and /3, hence actually 

«(1) = /«(2) = 0. 

Multiplying (33) by a^, we have 

3 

Integrating one more time, G is obtained as 

G = {G^i)Us + (^(2)) a'^+(G'(3)M3 + G'(4)) a^+(G'(5)M3 + G^q)) a^+(G'(7)M3+G'(8))a+G'(9). 

(35) 

The coefficients are found explicitly from -ft'292) but we omit these expressions. 

It follows that for P 7^ integrable evolution equations with nontrivial p(3) are 

of the form 

5 /5DF \ 

ut = a^u^ + - a^a'^ul + i-—^a^ + 5 {Da - 03^4) a" JM4 + (G'(i)M3 + ^(2)) a^ 

11 



+ (^(3)^3 + (^(4)) a^+ (^(5)^3 + ^(6)) a^+ {G(7)U3 + G(^s)) + ^(9), (36) 

where a and P are given respectively by (21) and (23), the G'(i), i = 1, . . . , 6 
are certain functions of a, (3 and 7, G^r) = G['jo)P-, where ^(70) is a function of 
X, t, u and Ui and G'(8), ^(g) are functions of x, t, u, Ui and U2. 

Proposition 3.1. Let Ut = a^u^ + i? be an evolution equation where a and B 
depend on x, t, u, Ui, i = 1, . . . , 4. If the canonical densities p(j), i = —1, . . . , 5 

are conserved, and q^^-' and „J|^ are nonzero, then integrable evolution 
equations are of the form in (36), where a and P are given respectively by (21) 
and (23), the G(j), i = 1,...,6 are certain functions of a, (3 and 7, ^(7) = 
G(jo)P, where G(jo) is a function of x, t, u and Ui and G(8), G(9) are functions 
of X, t, M, Ml and ^2- 

3.2.4 A new solution 

The equation (36) obtained above gives the form of the candidates for inte- 
grable fifth order quasi-linear equations. In this section we assume that a is 
independent of x, u and ui and obtain an equation for which all of the canonical 
densities p(j), i = — 1, . . . , 5 are conserved. 

Proposition 3.2. Let Ut = a^u^ + i? be an evolution equation where a and 
B depend on -Uj, i = 2, . . . , 4. If the canonical densities p(i), i = —1, . . . , 5 are 
conserved, and ^^^^^ and „ ^f' are nonzero, then 

Ut = a^U5 + Bul + Gu4+ [G(i)a^ + G(^s)0-^ + G(5)a^ + G(7)a)U3 + G(^g), (37) 

where 

P + /32 



a 



{aul + f3u3 + -f)-^/\ a = a^o)P\ /? = 0, 7 



4a 

P = {fiul + Z/-U2 + /t)""*^, g = AjlK — Z/^, 



5 / 15 5 \ 

B = -- a(o)P^a'^M3, "^ = ( ~ Y '^(0)'^^ + o^j '^^'""^ + '^) (^^us, 

45 27 / 40 

^^'^ = 128 "(0) ^"' ^^^ " ^^^ ' ^(') = 16 ""(0) ^^' ( - ^^ + Y^. 



/ lo \ 

G(5) = «(o) -P""^ ( g-P — ^/ij , ^(7) = G(jo)P, G(^Q) = G 



and /i, z/, K, a(o), G(7o) and G'(9o) are constants. 

It has been checked that all even canonical densities are trivial and all 
odd canonical densities up to p(5) are conserved. Thus the equation above is 
integrable in the sense of admitting a formal symmetry [1] . We have observed 
that all the even canonical densities up to p(4) are trivial while the odd canonical 
densities are nontrivial. 
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4 Discussion of the Results 

In the course of this study, the existence of non-polynomial integrable equations 
came as surprise, because we were expecting that the polynomiality in top 
three derivatives would be generalized to lower orders. We then remarked that 
this class seems to be related to the essentially nonlinear class of third order 
equations [1], 

Ut = H= [Aiul + A2M3 + A^y^'^{2A^u^ + A2) + A4, 

where Ai = Ai{x,u,Ui,U2). This similarity is suggested by the fact that the 

conserved density P(-i) is the same for both equations. Note that a scaling of 

the Aj's result in a scaling of H; replacing A^ by (2/P)^/^Aj, i = 1,2,3 one can 

obtain 

dH , 

Hence P(_i) is the same for the essentially nonlinear third order class and the 
fifth order quasi-linear non-constant separant class. 

Then, in order to see whether these equations were part of a hierarchy, 
we studied the top level terms of level homogeneous odd order equations up 
to order 13; surprisingly we have seen that at each order, these flows were 
polynomial in Uk, k > 3 oi the form Ut = a^Um + . . . , i.e, with separant a"^. 

Finally we note that the classification of essentially nonlinear class of equa- 
tions is discussed in [8] where it is suggested to use a formulation where the de- 
pendent variable is replaced by the canonical density P(-i). With this method, 
the classification of the essentially nonlinear equations is partly completed in 
[11]; in particular the classification of the equations with Ai = corresponding 
to Q;(o) = in our case is obtained, while the case Ai 7^ is still open. 

Acknowledgment 

This work is based on the Ph.D. thesis of the first author and it is partially 
supported by the Turkish Scientific and Technical Research Council. We thank 
Professors Metin Giirses and Jing Ping Wang for valuable discussions; for their 
suggestions to concentrate on specific solution (MG) and for pointing our recent 
work by Heredero (JPW). 



13 



Appendix A. Canonical Densities 

P(-i) = a'\ (A.l) 

P(o) = F4/F5, {A.2) 

2 1 

P(i) = 2 (D^a) - (Da)2a-i + 2 (Da) a-^F4 - - {DF4) a'^ + - a'^F^^ 

5 5 

-Aa-^F4^ (A3) 



P(2) = -2 (D^a) 0-4^4 + 12 {Dafa-^Fi - 4 (L)a) (£1^4) a"^ 

+3 (Da) a^^Fs - — {Da) a'^F^^ - \ {DF3) a'^ + \ (D^F^) a~^ 



+ g {DF,) a-'F, + I a-'F, - ^ a'^F^F, + -^ a~''F^, (AA) 



P(3) = --{D'^a)a^ - — {D^a){Da)a-{D^a)a-^F^ 
5 5 

12 19 

(D^afa + — (D^a) (Da)^ + 2 (D^^) (/)a) 0-4^4 

5 5 

-- (D^a) (DF,) a-3 - - (D^a) a~^Fs + - (D^a) a'^F^ 

5 5 25 

3 14 

-- (Da)^a-i + 2 {Dafa-^^F^ + — (Da)^ (DF4) a^^ 

5 5 

+ 1! (/)a)2 a-^F3 - — {Daf a~'F,^ - — (Da) (DF,) a'^ 
5 25 5 

— (Da) {D'^Fi) a~3 + {Da) (DF^) a~^F^ + 3 {Da) a'^F2 

5 25 

21 49 3 

_ (^Da) a-^FsF, + — (Da) a'^^F^"" - ^ {DF^) a'^ 

5 25 5 

+ \ (D'Fs) a-' + A ^DFs) a-'F, + i (D'F,) a'' 
5 25 5 

-1 (^2^4) a-^F4 - ^ (DF,)' a-^ + | (DF4) a'^F, 

--4 (^i"4) 0-^2^4' - I a-V(_i) + I a-^Fi - 4 a^'F^F^ 
125 5 5 25 

_1 a-^i^gS + ^ a-'^FsF^ - —a~'''F^\ {A.5) 

25 125 625 ^ ^ 
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4 fifi Ifi 

p,^. = - (D*a) a-^F^ - — (D^'a) (Da) a~^F^ + — (D^a) (DFA a'^ 

5 5 5 

-iD^a) a'^Fs + - (D'a) a-'F,' - — (D^a? a'^F^ 
5 5 

+ ^ (D^a) {Dafa-^F^ - — {D^a) (Da) (DF^) a'^ 
5 5 

-8 (D^a) (Da) a'^F^^ - - (D^a) (DF3) 0"^ + — (D^a) (D^FA a"' 

5 5 

+ — {D\) (DF^) a-'F^ - - (D^a) a'^F^ + — {D^a)a-'F^F^ 
25 5 25 

44 264 1 Q2 

_il m^a) a-'^F^^ - — iDa)^a-''F^ + — iDa)\DFAa-^ 
125 5 5 

+3 (Da)3 a^^Fg + — (Da)^ 0-9^4" + - iDa)"" {DF^)a-'' 
5 5 

-12 (Da)' (^2^4) a-3 - — {Da)^ (DF^) a'^F^ + — {Daf a'^F^ 

-^ (Da)' a-'F,F, + ^ (Da)' a-i3F4^ - ^ (Da) (DF,) a'' 

-^ (Da) (D^Fs) a-' + ^ (Da) (DF3) a-'F4 + 2 (Da) (D='F4) a"' 
5 25 

64 44 

+ — (Da) (D^F^) a-'F^ + — (Da) {DF^f a~' 

+ — {Da) (DF^) a-'Fs - -— (Da) {DF4) a~^^F^^ + 2 {Da) a~^Fi 

12 6 66 

-— {Da) a-'F^Fi - - {Da) a"' F^^ + — {Da) a-'^F^F^^ 

ZiO 

'''' {Da) a-^^Fi^ - - (DFi) a~^ + — (DF2) a-^F4 



125 ' ' 5 ' ' 25 

-i (D^Fs) a-i - ^ (DF3) (DF4) a-6 + ^ (DF3) a-^Fs 

^^ (DF3) a-"F4' - — (D^Fi) a-i - — (D^Fi) a-^^F^ 



125 ' ' 25 ' ' 25 

-| (D'F4) (DF4) a-6 - I (D'F4) a-«F3 + ^ (D'F4) a-i^F4' 

24 8 48 

+ 125 (^^4)'a""F4 + -(DF4)a-6F2 - — {DF^)a-^^F^F, 

+ ^{DF^)a-''F4^'' + 2(Da(_i)) + - a'^Fo + — a' V(o) 



P(5) = 4(a V(i) -p(i)C7(_i)), (A.7) 

where Fm = dF/dum, m = 1, . . . , 5 and <J(i),i = —1, 0, 1 are differential poly- 
nomials such that (p(j))^ = Da(^iy 
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